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Abstract. This series of papers is dedicated to a new technique to select galaxies that can act as standard rods and standard 
candles in order to perform geometrical tests on large samples of high redshift galaxies to constrain different cosmological 
parameter The goals of this paper are (1) to compare different rotation indicators in order to understand the relation between 
rotation velocities extracted from observations of the HQ■/^6563A line and the [OII]/i3727A line, and (2) to determine the scaling 
relations between physical size, surface brightness and magnitude of galaxies and their rotation velocity using the SFI++, a 
large catalog of nearby galaxies observed at I-band. A good correlation is observed between the rotation curve-derived velocities 
of the Ha and [Oil] observations, as well as between those calculated from velocity histograms, justifying the direct comparison 
of velocities measured from Ho- rotation curves in nearby galaxies and from [Oil] line widths at higher redshifts. To provide 
calibration for the geometrical tests, we give expressions for the different scaling relations between properties of galaxies 
(size, surface brightness, magnitude) and their rotation speeds. Apart from the Tully-Fisher relation, we derive the size-rotation 
velocity and surface brightness-rotation velocity relations with unprecedentedly small scatters. We show how the best size- 
rotation velocity relation is derived when size is estimated not from disc scale lengths but from the isophotal diameter (-23.5, 
once these have been corrected for inclination and extinction effects. 



1. INTRODUCTION 

With the new generation of telescopes and instruments, such 
as multi-object spectrograph and large field of view cameras, 
it has become possible to perform large, deep redshift surveys. 
Using these data, the structure and geometry of the Universe 
can be studied through new channels, namely geometrical tests 
such as the angular diameter test (angular size - redshift re- 
lation), the Hubble diagram (magnitude - redshift relation) and 
the Hubble test (count - redshift relation). These tests are gener- 
ally difficult, requiring differentiation between the natural evo- 
lution of galaxies and the effects of geometry. Large redshift 
surveys ease the task in that respect. 

To perform the geometrical tests mentioned above, a pop- 
ulation of objects that can be tracked through redshift needs 
to be defined. More specifically, the angular diameter test re- 
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quires a standard rod to be identified, the angular size of 
which is then being measured over the redshift range of in- 
terest. The objects taken to serve that pur pose could be any- 
thing fr om galaxies (e .g. Sandagj, 1972 ). clusters of galax- 
ies (e.g. IS ereno, '2003t ICooravl 19981) . or dark matter halos 
dCoorav et a l.. 200 li). Regardless of its nature, a standard rod 
should satisfy simple criteria: its size should be measured re- 
liably up to high redshifts, and it should be observable in the 
local Universe to provide calibration. In this study, the relation 
between the physical size of a spiral galaxy and its rotation 
velocity (or global profile width) is used. By virtue of this re- 
lation, selecting a population of objects with a given rotation 
velocity is equivalent to selecting them based on their physical 



In other papers of this series (*lV Iarinoni et al. , l2006allbl 



thereafter Paper I and Paper III), the angular diameter test is 
performe d on a dataset of prec ursor Vimos/VLT Deep Survey 
(VVDS) (le F evre et all bOOl *) galaxies, using the linewidth- 
diameter relation to select standard rods. The test is used to iso- 
late the effects of disk evolution with redshift, and to constrain 
the values of cosmological parameters. Though the VVDS pro- 
vides excellent data at high redshift, the angular diameter test 
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highlights one of its shortcomings: Hmited volume coverage at 
low redshift and therefore the lack of a comparison sample in 
the local Universe. 

In this context, the goal of this paper is tw ofold. First, us 



inga large sample of nearby galaxies (SFI++: ISpringob et al. 



20071) the physical properties of the standard rods and standard 



candles used to perform the geometrical tests are established, 
free from any evolutionary biases. This is especially useful to 
calibrating the angular diameter relation and to separating the 
effects of galaxy evolution from those of geometrical varia- 
tions (see Paper I). By establishing clearly the scaling relations 
between the size, surface brightness and magnitude of these 
nearby galaxies observed at I-band and their rotation velocity, 
we establish the calibration that will allow us to select standard 
rods/candles simply by tracing through redshift galaxies with 
a given rotation speed. Having a strong handle on the struc- 
tural parameters of disc galaxies in the local universe gives the 
unique opportunity of tracing over time the evolution of these 
properties for disc galaxies hosted in dark matter halos of the 
same mass, if rotation velocity is used as a proxy for halo mass, 
and compare them with predictions made un der the hierarch i- 



cal scenario for the growth of structures (e.g. lMo et al.Lll998b 



The second goal of this paper is to compare the rotation ve- 
locity indicator used for the high redshift data to those used in 
the local Universe, cross-calibrating rotation indicators used at 
different redshifts. This last point is of relevance since rotation 
information for spiral galaxies can be obtained through the ob- 
servation of various spectral lines, the choice of which varies 
with the redshift of the sample. For galaxies with z < 0.1, the 
HI 21 cm line is an excellent candidate. The HQf/16563A line 
is also frequently used for galaxies with low to moderate red- 
shifts. However, the Ha line is quickly redshifted into the near- 
infrared and becomes unavailable to ground observers using 
optical telescoptes for galaxies with z > 0.4, even though with 
new near-infrared spectrographs such as SINFONI it is now 
possible to obtain Ha rotation c urves for high redshift galax- 
ies (e.g. IForster Schreiber et all i2006). Even with the advent 



of such instruments, most large studies of galaxies at high red- 
shift rely on the [OII]i3727 A hne, including th e VVDS and 



the DEEP2 Redshift Survey (iDavis et al.L l2001h . In order to 



compare sets of local and distant galaxies, it is therefore neces- 
sary to understand how rotation velocities extracted from these 
different hnes relate. It has already been shown that velocity 
widths derived from HI 21cm and Ha observations are in ex - 
cellent agreement (see for example Courteau, 1997; Vogt et al. 



[20041 The coiTelation between the [Oil] and HI Hne wi dths has 
also been investigated ( Kobulnicky & Gebhardt . 2000h . Using 
a sample of 22 nearby late-type spirals they find [Oil] widths 
to be accurate to within 10% for galaxies with a roation width 
of 200 km s ', which is comparable to the overall scatter in 
the local Tully-Fisher relation. However, they conclude that 
the uncertainties go up to about 50% for galaxies with widths 
< 150 km s Here, a sample of 32 spiral galaxies with 
0.155 < z < 0.25 is used to compare velocity widths from 
the Ha/I6563A and [OII]i3727 A Hnes (hereafter. Ha and 
rOIII). Even though these galax ies are more distant than in the 



a better estimation of the rotation velocities. Our data also al- 
lows for a direct comparison between Ha and [Oil], therefore 
bridging nicely the gap between low and high r edshift sam- 
ples. Contrary to the Kobulnicky & Gebhardtl ( 2000) study, the 
galaxy sample used here is restricted to late spiral types, there- 
fore we will not address the question of possible biasses due to 
varying morphologies. 

In ^a description of the galaxies observed and of the data 
reduction process is given, while in ^we describe how veloc- 
ity widths were extracted from the data and how they compare 
for both sets of emission lines. In 34.1l we presente a sample of 
nearby galaxies used as calibrators for the purpose of the an- 
gular diameter test (Paper III) and discussion of our results are 
given in ^ Calculations were done assuming Hq = 70 km s"' 
Mpc-\ Q.M = 0.3, = 0.7. 



2. DATA 

2.1. Target Selection and Observational Setup 

Using the Hale 5m telescope of Mount Palomar during the 
course of three short observing runs between 2003 March and 
2004 February (total of 12 nights), long-slit spectra were ob- 
tained for a sample of spiral galaxies. The galaxies were cho- 
sen in the area covered by t he Data Release 2 of the Sloan 
Digital Sky SurveJ^ (SDSS, IStoughton etall l2002b . and se- 
lected to have inclinations in excess of 50 degrees, Petrosian 
radii containing 90% of the light in r band larger than 4" to 
avoid AGN-dominated systems, and [OII]/l3727A equivalent 
widths in the range of 10-35 A. The galaxies were also required 
to be in the redshift range 0.155 < z < 0.25. Galaxies with 
close neighbours or disturbed morphologies were rejected to 

avoid including interacting systems in the sample. 

Using the Double Spectrograph dOke & Gunnill982 l). both 
Ha and [Oil] emission lines could be observed simultaneously. 
The main properties of the galaxies are listed in Table [T] The 
first two columns are the AGC (Arecibo General Catalog) and 
SDSS names, respectively (the J2000 coordinates of the objects 
can be directly extracted from their SDSS name). The third col- 
umn is the heliocentric systemic velocity of the galaxies, as 
measured from the Ha line. All other parameters (redshift, I- 
band magnitude, disc scale-length, inclination and position an- 
gle measured counterclockwise from North) are from the SDSS 
database. Unfortunately, about 50% of the observing time was 
lost to bad weather conditions, forcing the final sample pre- 
sented in this paper to contain only the 32 galaxies for which it 
was possible to get rotation curves in both Ha and [Oil] . 



Kobulnickv & Gebhardt (|200(^ sample, the scatter in the ro- 



tation curves data points is significantly smaller, resulting in 
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The red camera of the Double Spectrograph was used with 
a 1200 lines mm"' grating (7100A blaze) set at an angle of 
47° 18' in order to cover the wavelength range of 7520-8185 A. 
The spectral resolution is 0.65 A pixel ' with 24fim pixels and 
the spatial scale is of 0'.'468 pixel"'. The CCD camera had a 
size of 1024x1024 pixels, and a l"xl28" sHt was used (ex- 
cept for a few galaxies where the 2"xl28" slit was used due to 
very poor seeing conditions). The blue camera was set to cover 
the range 3675-5140 A by using a 1200 lines mm"' grating 
with a 4700A blaze at an angle of 35°52'. Using the 2788x512 
chip {ISfim pixels), a spatial sampling of 0'.'390 pixel"' was 
achieved, as well as a 0.55 A pixel"' spectral resolution. 

The exposure time was of 3600 s for most galaxies, though 
it was extended to 4800 s for a few sources when the flux in the 
[Oil] line was very low. Sky transparency and seeing conditions 
varied greatly over the course of the observing runs, but most 
observations were made between 1" and 2" seeing. In order 
to align the slit with the major axis of the galaxies, the ring 
angle of the spectrograph was set using the position angles of 
the galaxies given in the SDSS database. For each galaxy, the 
ring angle used is given in Table[T] On each night of observing, 
standard calibration images were obtained: dome flats, twilight 
flats (for the blue camera only), bias frames and darks. Given 
observing conditions, no attempt was made for any photometric 
calibration. 

2.2. Data Reduction 

The data from both cameras were reduced independently using 
tasks f r om th e Image Reduction and Analysis Facility (IRAF) 
dTodvl 119931) . Images from the blue camera had to be re- 



duced with particular caution due to the presence of many bad 
columns and other defects in the CCD. First, all the images 
were trimmed to remove overscans and low sensitivity regions 
on the edges. Then flat fields were constructed separately for 
each camera and for each night using the dome flats, as well as 
the twilight flats in the case of the blue images. After this, the 
images were also corrected for bad pixels using a mask created 
from twilight flats and the effects of cosmic rays removed. 

The cleaned images were wavelength calibrated and cor- 
rected for s-shape distortion. This was done in the dispersion 
direction by using observations of a reference star that was po- 
sitioned at five diff'erent places along the slit, and for the spa- 
tial direction by using either strong sky lines (red camera im- 
ages) or lamp spectra (blue camera images). The correction was 
in some cases important, especially in the dispersion direction 
where it was generally on the order of a few pixels. Individual 
exposures of each galaxy were then appropriately shifted and 
stacked. For these combined images, the background was fit- 
ted in order to remove the sky lines. The process allowed for 
optimal removal of the sky lines falling near or on the galaxy 
emission lines of interest. The continuum was finally subtracted 
from these background corrected frames. 



3. COMPARING AND CALIBRATING ROTATION 
VELOCITY INDICATORS 

Using these reduced data, it was possible to extract rotation 
velocity information for the galaxies in this Palomar sam- 
ple (hereafter referred to as the zO.2 sample) using both Ha 
and [Oil] emission lines and two different techniques: rotation 
curves and velocity histograms. We here describe these two dif- 
ferent techniques. 

3.1. Rotation Curve Extraction and Modeling 

A Gaussian was fitted along the dispersion axis at every pixel 
inside a box encompassing the emission line (either Ha or 
[Oil]). In the case of the [Oil] line where the signal-to-noise 
ratio (S/N) was often smaller than 3, Hanning smoothing was 
performed before the Gaussian fitting took place. Then, the ro- 
tation curves were calibrated, the measured wavelength con- 
verted to a velocity and the radial axis scaled to arcseconds. 
Following this, each point of the fit was examined and either 
rejected or included in the rotation curve. The criteria were a 
minimal S/N for the fitted Gaussian (4 for Ha and 3 for [Oil]) 
and a maximum width of 400 km s"'. The resulting rotation 
curves were then examined and some rejected points reintro- 
duced, while some others that satisfied the criteria were judged 
unrealistic and manually removed. These points for which a 
manual intervention was required were almost all either in the 
outer regions of the galaxy where the S/N drops significantly 
or contaminated by sky lines or cosmic rays that were not com- 
pletely successfully removed. Finally, the center of light of the 
galaxy was determined and used to center the rotation curve 
along the sp atial direction. Th is procedure is similar to the one 
presented in lVogt et al.l ( 120041) . 

The rotation curves were subsequently modeled by first 
folding about the kinematic center and fitting a Polyex model 
dGiovaneUi & Havnesl 120021) of the form: 



y(r) = fli(l-e"°^0(l+a3r). 



(1) 



Astronomy (MPIA), the Max-Planck-Institute for Astrophysics 
(MPA), New Mexico State University, University of Pittsburgh, 
Princeton University, the United States Naval Observatory, and the 
University of Washington. 



The parameter a i is related to the maximal velocity of the ro- 
tation curve and 02 controls the exponential rise of the rotation 
curve and the position of the "elbow", the point where a tran- 
sition occurs between the steeply rising inner part and the flat 
outer region. Finally, controls the slope of the outer part of 
the rotation curve. Due to the small radial extent of the emis- 
sion lines, the value of 03 is meaningless in most cases and has 
been manually set to zero. The Ha rotation curves are shown 
for all galaxies in Figure [T] and the [Oil] rotation curves in 
Figures |2l The solid lines are the Polyex fits, while the vertical 
dotted lines show the position of rgs, the radius containing 83% 
of the light of the galaxy. 

In Figure [3] the Polyex fits to both the Ha and the [Oil] 
rotations curves are shown together when sufficient data was 
available in both spectral lines to allow for the fitting process. 
In all the cases where the [Oil] line was bright enough for 
the rotation curve to be traced out to the region where it starts 
flattening out, the agreement between the two sets of rotation 
curves is excellent. In the other cases, the trend is for the [Oil] 
rotation curves to be slighdy shallower. 
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Fig. 1. The Ha rotation curves and their Polyex fits for the 32 spirals of the zO.2. The dotted lines show for each galaxy the 
position of rgs, the radius containing 83% of the light, and the solid line is the Polyex fit (when sufficient data was available). 
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Fig. 2. The [Oil] rotation curves and their Polyex fits for the galaxies which had detectable LOII]/l3727A emission. The sohd hne 
represents the Polyex fit (when sufficient data was available), and the dotted Unes show for each galaxy the position of rga, the 
radius containing 83% of the Ught. 
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Fig. 3. Polyex fits to the Ha (solid line) and to the [Oil] (dashed line) rotation curves presented in Figures[T]and|2] when both fits 
are available. Again, the vertical dotted lines show for each galaxy the position of rga. 
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AGC# 


Name 


Vo 


Zs DSS 


m, 


ro 


i 


PA 






kms" 1 




mag 


kpc 


dee 


deg 


101738 


SDSS J001107.7+003551.7 


46600 


0.155 


17.81 


4.24 


66 


8 


101701 


SDSS J003123. 8+141313. 8 


49444 


0.165 


17.15 


6.93 


67 


116 


112393 


SDSS J014451.5+133616.9 


46711 


0.156 


17.31 


6.18 


54 


111 


130796 


SDSS J031818. 8+010750.4 


48640 


0.162 


17.68 


3.29 


51 


87 


431415 


SDSS J032434.8-07 1050.6 


49673 


0.166 


17.32 


7.82 


64 


356 


130826 


SDSS J033331. 8+010717.1 


53337 


0.178 


17.48 


6.92 


61 


158 


431418 


SDSS J040000.2-053852.3 


50138 


0.167 


17.36 


7.03 


66 


51 


171385 


SDSS J073235.0+374734.4 


65413 


0.219 


16.94 


9.09 


62 


331 


171389 


SDSS J074608.6+395658.5 


47142 


0.157 


17.46 


5.54 


49 


343 


171394 


SDSS J075555.0+364057.7 


55351 


0.185 


17.01 


7.60 


49 


103 


181489 


SDSS J080438.6+463536.6 


60671 


0.202 


17.12 


6.00 


55 


12 


181487 


SDSS J081345.7+473552.5 


54478 


0.182 


17.40 


6.28 


61 


322 


181490 


SDSS J081750.4+48151 1.3 


55617 


0.186 


17.30 


7.82 


56 


35 


181488 


SDSS J084304.2-000737.6 


49300 


0.165 


17.40 


8.79 


63 


7 


191728 


SDSS J091957.0+013851.6 


52850 


0.176 


16.90 


5.51 


51 


20 


191729 


SDSS J092050.9+575758.5 


51558 


0.172 


16.94 


6.01 


46 


117 


191731 


SDSS J094412.5+020206.5 


47453 


0.158 


17.17 


6.30 


60 


99 


191732 


SDSS J100047.6+020406.6 


55778 


0.185 


16.71 


5.02 


61 


90 


501752 


SDSS J100757.3-005612.3 


46402 


0.155 


17.08 


4.38 


58 


179 


202041 


SDSS J104215.5+010508.3 


58452 


0.195 


17.46 


7.47 


60 


104 


501753 


SDSS Jl 101 17.7-010014.9 


47277 


0.158 


17.39 


7.15 


66 


36 


212861 


SDSS Jl 10936.9+010957.8 


63918 


0.213 


17.42 


9.08 


60 


95 


511373 


SDSS Jl 14001.2-004624.3 


56368 


0.188 


17.36 


5.25 


58 


107 


212862 


SDSS Jl 14919.0+001806.3 


48216 


0.161 


17.14 


6.40 


53 


102 


511374 


SDSS Jl 15640.0-005 128.6 


50585 


0.169 


17.13 


6.66 


55 


122 


521390 


SDSS J123824.5-013339.8 


56357 


0.189 


17.03 


6.06 


54 


41 


224090 


SDSS J124819.3+672312.2 


50790 


0.170 


17.38 


7.56 


51 


26 


224092 


SDSS J125522.8+672306.8 


49401 


0.165 


17.51 


6.73 


55 


173 


242064 


SDSS J141054.1+013714.3 


51706 


0.173 


17.32 


8.46 


49 


129 


242065 


SDSS J141322.7+644404.7 


54763 


0.183 


17.51 


5.04 


53 


79 


242067 


SDSS J142846.2-000245.6 


55232 


0.184 


17.23 


7.75 


58 


242 


332257 


SDSS J234325.9+005605.2 


67105 


0.224 


18.20 


7.45 


64 


51 



Using these fits, the rotational velocity of each galaxy 
was determined, both for its Ho- and [Oil] lines. The veloc- 
ity adopted for a given rotation curve is the value of the Polyex 
fit at a radius corresponding to rgj. The values for r^T, are taken 
to be 3.5 X ro, where ro is the disk scale length. Since all the 
galaxies in the sample were taken from the SDSS, the values of 
physical parameters were taken from their database. The scale 
lengths are from an exponential fit performed on the i-band im- 
age of each galaxy. These values are presented in column 5 of 
Table [1] The value of for each galaxy is also represented in 
Figures[T]and|2]by the vertical dotted lines. Also, Tables|2]and 
[3]present the values of the Polyex coefficients and the rotation 
velocity at r83, Vrot{r%-i), for the Ha and [Oil] data, respectively. 

In Figure |4] we compare the rotation velocities that come 
from t he Hg and [OH] r otation curves. Following the prescrip- 
tion of llsobe et al.l (11990 ). a best fit to the data was obtained 
by linear regression using an ordinary least-squares bisector 
method (also called "double regression"). The two outlying 
data points were ignored during the fitting. The solid line in 
Figure |4]is the best-fit relation and the two dashed lines are the 
minimum and maximum slopes allowed by the best-fit param- 
eters. The dotted Une is the 1 : 1 relation. The standard deviation 



of the residuals of the best-fit (shown in the bottom panel) is 
28.6 km s"'. 



3.2. Velocity Histograms 

In order to cross-correlate the results with techniques used at 
higher redshift, measurements of the rotation velocities using 
velocity histograms were also made. These velocity histograms 
were built by collapsing the two-dimensional images along the 
spatial direction in order to form one-dimensional spectra. To 
determine the region of the image that will be collapsed, the 
central position of the continuum along the spatial direction as 
well as its width were determined by fitting a simple Gaussian. 
This region needs to be large enough to include all the emis- 
sion from the galaxy, but not unnecessarily wide to avoid in- 
creasing the noise in the one-dimensional spectra. Since some 
of the galaxies have their emission lines extending slightly be- 
yond the continuum, everything within 5cr of the center of the 
continuum was selected and added together, thus producing an 
intensity weighted spectrum. 

A Gaussian was fitted to each of the Ha histograms. The 
full width at half maximum (FWHM) of that Gaussian was 
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Table 2. Galactic Rotation Properties, from the Ha Emission Line 
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converted into an estimate of the rotation velocity of the galaxy, 
(Trot, via: 

(FWHM^ - FWHM2)1/2 ^ 1 

~ 2.355 Aq{\ +z) ' sin(0 

where FWHM/ is the FWHM of the emission line histogram, 
FWHM, the instrumental broadening, Aq the rest-frame wave- 
length of the emission line and / the inclination of the disk. 
The inclinations were determined from the axial ratio of the 
disks as derived from the same SDSS fits that produced the disk 
scale lengths. In the case of the Ha lines, FWHM,- was obtained 
by fitting a Gaussian to a bright night sky line that falls near 
the emission line. The velocities obtained using this technique 
((Trot) are plotted against the velocities derived from the rota- 
tion curves (vrot) for the Ha emission on the top panel of Figure 
|5]and are presented in column 6 of Table[2] In Figure|5] we plot 
only the galaxies for which the rotation curve extends radially 
far enough to reliably determine a rotation velocity. The stan- 
dard deviation of the residuals of the o-roi(Ha) - v,oiiHa) from 
the 1:1 relation is 20.8 km s which represents an eiTor of 
~ 10% for a galaxy with a rotation velocity of 200 km s"'. 

For the [Oil] emission, one-dimensional spectra were pro- 
duced in the same manner, but this time cr„, cannot be calcu- 



lated by a single Gaussian fit since the doublet nature of the 
line needs to be taken into account. Instead, deblending of the 
two lines was done using the "splot" procedure in IRAK The 
expected positions of the two lines of the doublet were given 
as input, and the relative position of the two lines is maintained 
fixed during the fitting process. We farther required that the 
width of the two lines be the same (as it should be for phys- 
ical reasons), and that the amplitude of the two lines be in a 
reasonable ratio (not more than a fac tor of three between the 
two, which is similar to the criteria of ^ Kobulnickv & Gebhardt 
(2000)). Since the process relies significantly on subjective 
judgement on the part of the data processor, especially with 
regard to the size of the fitting region around the emission line, 
it was performed independently by two of us (AS, CM). The re- 
sults were compared and a quality flag was put on the values. A 
width estimate was judged reliable when the two values agreed 
within the spectral resolution of the observations (i.e. 0.55A). 
The widths were transformed into rotation velocities using eq. 
|2l This time, however, FWHM, is calculated using exposures 
of a FeAr hollow cathode calibration lamp, since there are very 
few bright sky lines over the relevant wavelength range. The 
calculated velocities (cr,.,,,) are plotted against the rotation curve 
derived velocities for Ha on the bottom panel of Figure |5] and 
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given in the last column of Table [3] The rms of the residuals of 
the crrot({OII}) - Vrot(Ha) relation is 52.7 km s"' , which at any 
rotation width represents an error larger by a factor of ~ 2.5 
than the relation for Ha alone. 

3.3. Definition of the Rotation Classes 

When comes the time to perform the geometrical tests, the 
VVDS rotation velocities are used to separate the galaxies into 
two samples, one of slow and one of fast rotators. This is done 
in order to trace out the evolution of two different populations 
of galaxies by using the angular diameter test (see Paper I and 
Paper III). These populations are defined as the "slow rotators", 
which have < v,o, < 100 km s"' and the "fast rotators" with 
100 < Vroi < 200 km s"'. The analysis of the rotation curves 
and of the velocity histograms for the zO.2 sample has shown 
that even considering the scatter in the relation between the dif- 
ferent velocity indicators and the different spectral lines, simi- 
lar criteria on rotation velocities can be applied across the dif- 
ferent data sets to define the populations of objects that can act 
as standard rods or candles to perform the geometrical tests un- 
der the approach suggested in Paper I. Because of the restricted 
size of the zO.2 sample and the very small number of "slow ro- 



tators" with reliable velocity measurements, these galaxies will 
not be used directly to perform the geometrical tests in Paper 
III. However, given the bridge that we have just established 
between low and high redshift velocity indicators, one could 
eventually use large data sets of nearby galaxies, such as SDSS, 
to complete the gap in redshift coverage between the very local 
SFI-H-i- catalog (see next section), and the VVDS data. 



4. CALIBRATING THE STANDARD RODS AND 
STANDARD CANDLES 

In order to perform some geometrical tests such as the angular 
diameter test, not only will we need the well-defined samples of 
galaxies at 0. 15 < z < 1 .4 but also reference values for the local 
Universe (z ~ 0). Since the contributions of galaxy evolution 
and geometric signature are superposed, it is essential to de- 
termine the sizes and luminosities of the standard rods/candles 
for a local sample that is free from any evolution, and also the 
scatter in these relations. 
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Fig. 4. Comparison of the rotational velocities derived from the 
rotation curves for Hot and [Oil]. Top Panel: the solid line is 
the best fit obtained by an ordinary least-squares regression, 
excluding the two outlying points. The dashed lines are the 
maximum and minimum slopes determined by the errors on 
the parameters of the best-fit linear relation and the dotted line 
corresponds to a 1:1 relation. Bottom Panel: the residuals of the 
best-fit model. In both panels, the filled symbols represent the 
15 galaxies with the most secure rotation velocity estimates, 
while the open symbols are those galaxies with poorer rotation 
curve extent, either in Ha or [Oil]. 



4. t . /A Sample of Galaxies in the Local Universe 

The SFI-i-i- dSpringob et all l2007h . a sample of ~ 4200 lo- 
cal galaxies, is used for this purpose. They range in distance 
from to 300 Mpc, with most objects within 100 Mpc. SF14-+ 
builds on the earlier SFl (Spiral Field 1-band) and SCI (Spiral 
Cluster I-b and) samples published i n a series of papers in the 
1990s (e.g. ' Giovanelli et all ( Il997b : iHavnes et all d 19991) ). but 
includes a significant amount of new data, plus reprocessing of 
the old values. Since SFI-i-H- is a very large and homogeneous 
sample of nearby spiral galaxies, it is perfectly suited for the 
study of scaling relatio ns in disc galax i es. A n 1-band TuUy- 
Fisher template relation ( Tullv & Fisheil 1977h is pres ented in 
detail for the SCI subsample in'Oiovanelli et al.' ('199?'), and for 
the larger SFI+-I- cluster sample in Masters et al. (2006, there- 
after M06). Here, the values of half-light radius, absolute mag- 
nitude and surface brightness are used as calibration for high-z 
data. 



4.2. Disc Scaling Relations 

Under the hierarchical scenario for the growth of structure 
scaling relations between these three quantities and th e rota 
tion velocity of dark matter systems is predicted (iMo et al 



Fig. 5. Comparison of the rotational velocities derived from the 
rotation curves (vror) and from the one-dimensional collapsed 
spectra (cr,.^,) for Ha (top panel) and [Oil] (bottom panel). In 
each pannel the dotted line corresponds to a 1:1 relation, and 
the continuous lines a deviation of 15% from it. 



1998h . The most commonly used of these, which is also the one 
with the less scatter, is the magnitude-rotation velocity (Tully- 
Fisher) relation . 

4.2.1. The Magnitude-Rotation Velocity Relation 

The TuUy-Fisher relation for the SFI+-I- sample is shown in 
Figure |6] The best fit to the Tully-Fisher relation at I-band has 
been derived by M06 for a late-type subset of these galaxies, for 
which accurate distances can be obtained through membership 
in nearby clusters of galaxies. We will from now on refer to this 
subset of galaxies as the "template sample". The Tully-Fisher 
relation from M06 is: 

M-Slogh = -20.85 ±0.03 -(7.85 +0.1)(logv„,-Hlog 2 -2.5)(3) 

We adopt this expression as our magnitude-velocity relation. 
M06 also give a precise estimate of the scatter around that rela- 
tion. We recompute it here however for the complete SFI-n- 
sample, to account for the broader morphological types in- 
cluded and the additional scatter coming from generally larger 
distance uncertainties. We measure the scatter to be cr = 0.56 
mag over the entire range of rotation velocities. 

4.2.2. The Size-Rotation Velocity Relation 

Since one of the main goals of this series of papers is to perform 
the angular diameter test on a sample of high redshift galaxies 
by using the relation between the physical size of a disc and its 
rotation velocity to select standard rods, it is most important to 
derive a reliable expression for the z = size-rotation velocity 
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Fig. 6. Magnitude-rotation velocity relation for the sample of 
~ 4200 SFI++ galaxies. Absolute magnitude at I-band, cor- 
rected for extinction and corrected to face-on (^ Springob et al. 



20071) . plotted against rotation velocity. The solid line in each 
panel is the template Tully-Fisher relation derived by M06, us- 
ing only the subset of galaxies with the most accurate distances 
(top panel). The bottom panel includes all the SFl-i-H- galaxies. 
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Fig. 8. 



. relation for the template sample (top panel) and 



for the full SFI-i-H- sample (bottom panel). The quantity plotted 
is ?"23.5, the radius corresponding to the fi = 23.5 isophote, after 
corrections for inclination, foreground extinction and cosmo- 
logical stretch. In each case the solid line is the best bisector 
linear fit to the relation, the parameters of which are given in 
Tableg] 




1.8 2.0 2.2 2.4 

loe(v„i [km/sD 

Fig. 1. r - Vrot relation for the template sample (top panel) and 
for the full SFI-n- sample (bottom panel). The quantity plotted 
is r83, the radius encompassing 83% of the light. No attempt is 
made to correct rgs for inclination. In each case the solid line 
is the best bisector linear fit to the relation, the parameters of 
which are given in Table HI 



('" - Vrot) relation. ISpringob et al.l ( 120071) published two sets of 
sizes for the galaxies in the SFI-i-h- sample: rgB, the radius con- 
taining 83% of the light, and r23.5 the radius corresponding to 
the = 23.5 isophote. While values of the disc scale length rd 
are available for SFI-i-h galaxies, they tend to be unrehable and 



produ ce scaling relations with large scatter dGiovanelli et al 
19941) . In Figures |2] and [8] we present the r - Vrot for rg3 and 



'"23.5, respectively. The scatter around the best-fit bivariate lin- 
ear fit is smallest when r23.5 is used. The Pearson correlation 
coefficient of the relation is 0.92 for r23.5 and 0.67 for r%T,. We 



will therefore not consider r83 any further, but th e parameters 
of the best direct, inverse and bivariate linear fits dlsobe et al 



1990l) to the r83 



„ relation are given in Table |4l as well as 
those for r23.5 - Vroi- 

The smaller scatter in the r23.5 - v„,f relation is explained 
by the fact that r23.5 as opposed to r^j, is measured where the 
disc is close to optically thin, so it is possible to correct it for 
inclination effects. While no inclination correction has been at- 
tempted for r83, the values of r23.5 plotted in Figure[8]have been 
corrected for inclination to their face-on values using Equation 
7 of Giovanelli et al. (1995) . where a complete analysis of the 
inclination dependence of isophotal radii was done. We have 
further corrected these radii for foreground galactic extinction 
as: 



'23.5 



(4) 

1.086 

where ^ is the observed value coiTected to face-on. A,- the 
galacti c extinction at I-band derived from the DIRBE dust 
maps (ISchlegel et al.l Il998h . and the disc scale length. 



Both the corrections for inclination and extinction require a 
measurement of rj, but this quantity was not published by 



Springob et al.l ( l2007h for the SFI-n- sample, because less reli- 



able than the isophotal measurements of radii. We use these un- 
published values of to derive as a function of magnitude the 
ratio r23.5/rd, which is then used to scale the reliably-measured 
value s of 03.5 down to (this process follows the study of 



.Giovanelli et al.l (Il995h ). The small scatter in the r23.5 - Vrot 
can also be partly explained by the care that has been taken 
in deriving distances for the SFI+-1- galaxies. These distances 
were calculated by taking into account the peculiar velocity of 
each galaxy. For field galaxies, the peculiar velocities are de- 
rived from the systemic velocity of the galaxies and the abso- 
lute magnitude compared to the values predicted by the Tully- 



12 



Saintonge et al.: Geometrical Tests of Cosmological Models 11. 



a. 




1.8 3.0 2.2 2.4 

loe(v„u Tkm/sl) 



Fig. 9. 



relation for the template sample (top panel) and 



for the full SF1++ sample (bottom panel). The quantity plotted 
is rso, which has been computed as r^o - 1.68rrf, where the 
disc scale length was estimated by scaling the value of r23.5 by 
a factor depending on the magnitude of the galaxies. In each 
case the solid line is the best bisector linear fit to the relation, 
the parameters of which are given in Table ID 




Fig. 10. Surface brightness-rotation velocity relation for the 
sample of ~ 4200 SF1++ galaxies. Surface brightness at 1- 
band within rso, the radius containing 50% of the light, plotted 
against rotation velocity. The top panel uses only the smaller 
sample used by M06 to derive the template TuUy-Fisher rela- 
tion, while the bottom panel includes all the SFl-H-i- galaxies. 

plotted is the average surface brightness within r^Q, calculated 



Fisher relation. For galaxies in groups, the peculiar velocity 
is the wei ghted average of the p eculiar velocities of the group 
members dSpringob et al.L 120071) . The distances are further cor- 



M50 = niio, + 5 log(r5o) H- 2.5 log(2;r). 



(5) 



rected for the Malmquist bias. 

To calibrate the geometrical tests performed in Paper 111, 
the relation between the half-light radius r^o and rotation ve- 
locity needs to be derived. While it would be straightforward to 
calculate rso from rgs, we rather choose to calculate it from r23.5 
since the ras.s - Vrot relation is significantly tighter. Assuming 
that the discs obey an exponential profile, the half-light radius 
is simply r^o - LGBirj. Instead of using the directly measured 
but inaccurate values of r^, we use the values obtained by scal- 
ing down the inclination- and extinction-corrected values of 
r23.5 as explained above. The r^o - Vr„t relation is presented in 
Figure |9] and the parameters of the best-fit linear relations are 
in Table |4] This r^o - Vrot relation is central to the geometrical 
tests of cosmological models performed in Paper III. 

4.2.3. The Surface Brightness-Rotation Velocity 
Relation 

The surface brightness-rotation velocity (ju - Vrot) relation is 
the least often used of these three scaling relations, as it in- 
trinsically has more scatter than both the Tully-Fisher and the 
r - Vrot relations. This larger scatter is due to a stronger depen- 
dency of the ^ - Vrot relation on the spin parameter of the disc /I; 
while the M - Vrot relation is independent of A, this parameter 
appears linearly in the r - Vrot re lation and quadratically in the 
relation dMo etaUll998i) . 
In Figure[TO]we show the yU - Vrot relation both for the tem- 
plate sample and for the entire SFI-n- relation. The quantity 



The best fit relation to the entire SFI-i-H- sample, using an ordi- 
nary least squares bisector method, is: 

A(50 = (29.67 + 0.22) - (4.43 + 0.09) log(v™,) (6) 
The scatter of the residuals of this best fit relation is cr = 0.78. 

5. CONCLUSIONS 

The main goal of this paper was to provide calibration in two 
different ways for the angular diameter test that we performed 
on data from the VVDS: (1) by directly comparing velocity 
indicators commonly used in the local and distant Universe, 
and (2) by using a large sample of nearby galaxies to provide 
zero points for the radius - rotation velocity, magnitude - ro- 
tation velocity and surface brightness - rotation velocity re- 
lations. Previous studies have focused on comparing the Ha 
or [Oil] rotation velocities with HI data for nearby galax- 
ies. Kobulnicky & Gebhardt (2000) have used a sample of 22 
nearby galaxies to compare [Oil] and HI linewidths. Their 
galaxy sample was smaller and less homogeneous that the one 
used in this study. They found the two sets of velocity esti- 
mates to be consistent within 10% overall, with the [Oil] un- 
derestimating the HI rotation velocity by up to 50% in a few 
extreme cases. In the local Universe, the consistency of the HI 
an d Hg rotation velocities was also e stablished (see for exam- 
ple 'C ourteaul 1 19971 l\bgt et al.ll200 4'). Most recently, 'Catinelia 
(20051) has shown that on average HI widths are larger than Ho; 
ones by ~ 10% for a galaxy with w = 100 km s"' and that op- 
tical rotational widths measured from velocity histograms are 
affected by systematic biases and therefore are less reliable that 
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Table 4. Coefficieiits of best direct Jiiverse aiidbivariate Imear fits to the r-vvo/ relatiori, log(r[kpc h ']) = a+felog(vro,[km s ']). 



Relation 


Cldij- 


bdir 


a,,,,, 


bill 




ai,i 


bbi 


o-6,(log(r[kpc h"'])) 


Template galaxies 


















83 ~ ^Vo/ 


-0.36 ± 0.08 


0.56 ± 0.04 


-3.56 ±0.29 


2.04 ± 


0.13 


-1.43 ±0.06 


1.06 ± 0.03 


0.16 


'■23.5 - Vrot 


-1.35 ±0.06 


1.03 ±0.03 


-2.27 ± 0.07 


1.45 ± 


0.03 


-1.77 ±0.06 


1.22 ± 0.03 


0.10 


rso - v„, 


-1.06 ±0.05 


0.73 ± 0.02 


-2.16 ±0.09 


1.24 ± 


0.04 


-1.54 ±0.04 


0.96 ± 0.02 


0.11 


All SFI++ 


















''83 - VVo/ 


-0.78 ± 0.03 


1.05 ±0.01 


-2.12 ±0.02 


1.06 ± 


0.01 


-1.45 ±0.07 


1.05 ± 0.03 


0.14 


''23.5 - 'Vrot 


-1.64 ±0.02 


1.15 ±0.01 


-2.06 ± 0.03 


1.34 ± 


0.01 


-1.84 ±0.02 


1.24 ±0.01 


0.07 


rsO - Vro, 


-1.39 ±0.02 


0.88 ± O.OI 


-2.09 ± 0.04 


1.20 ± 


0.02 


-1.71 ±0.02 


1.02 ±0.01 


0.08 



those derived by using the full spatial information through the 
fitting of the rotation curves. 

For this study, we first analysed the rotation curves for all 
the galaxies in our sample and concluded that the rotation ve- 
locities derived from the [OII]/13727 A line translate directly 
into the velocities obtained from Ha rotation curves, given that 
the sampling of the [Oil] rotation curve extends far enough to 
reach the region of constant (or slowly increasing) velocity. 
Even under this condition and as shown in Figure HI there is 
scatter in the relation of about 30 km s"'. The next step was 
to compare the techniques used to determine rotation velocities 
in the local Universe (Ho- rotation curves) and at high redshift 
([Oil] velocity histograms). Once again, there is a direct rela- 
tion between the two but significant scatter on the order of 50 
km s"' (see Figure|5]). 

We conclude that all the rotation velocity indicators studied 
here give comparable results, but before combining samples 
analysed with different methods one should be aware of the 
important scatter in the different relations. This scatter is due 
mostly to the poorer spatial extent of the [Oil] rotation profiles 
compared to the Ha ones obtained with the same integration 
time, and can be as large as 25% for a galaxy with w - 200 km 
s"' or 100% if w = 50 km s"' in the case of the cr,o,([OII])- 
v,of(HQ') relation (Fig|5]l. Even considering this large scatter, 
the fact that there is no systematic deviation from a 1 : 1 relation 
allows us to use both methods in their respective context to 
divide galaxies in classes based on rotation velocity, as required 
to perform the geometrical tests with the high redshift sample 
(see Paper III). In Paper I it was showed that this is a viable 
way of selecting standard rods for the purpose of the angular 
diameter test. Since one cannot perform the test only with high 
redshift data, it was of great importance to establish that a set 
of local calibrators can be reliably compared to it, even if the 
rotation measure is different. 

Once it is established that the various rotation indicators 
can be reliably compared, the next critical step is to derive the 
scaling relations between luminosity, size, surface brightness 
and rotation velocity for a set of nearby galaxies, free from any 
evolution. We used the SFI+-I- catalog, a very large homoge- 
neous sample of nearby spiral galaxies with I-band photometry 
and rotation information obtained either through optical long- 
slit spectroscopy or HI measurements. We have particularly fo- 
cussed our attention on the r - Vror relation which is central 
to the cosmological tests performed in the other papers of this 
series. We have shown that by carefully correcting sizes for in- 
chnation and extinction effects as well as obtaining accurate 



distances to the galaxies, the r23.5 - Vr„t relation can be almost 
as tight as the Tully-Fisher relation, as judged by their Pearson 
correlation coefficients. The relation derived using r23.5 as the 
size measurement is considerably tighter than the ones obtained 
using the disc scale length rd or a radius defined to contain a 
given fraction of the light of the galaxy (in our case r^^ or rso). 
The next phase for this study on scaling relations will be to see 
what constraints our r - Vrot relation imposes on models for 
disc galaxy formation, for example the range of spin parame- 
ters its tight scatter allows. These scaling relations also offer ex- 
citing possibilities for the study of disc galaxy evolution. Since 
they provide accurate values for luminosity, size and surface 
brightness as a function of rotation velocity, measurements of 
the properties for galaxies with the same rotation velocity (i.e. 
hosted in a dark halo of the same mass) over a large redshift 
range can provide in a new and unique way information about 
the evolution of galaxies and a mean of testing predictions from 
models under the hierarchical scenario for the growth of struc- 
tures. A pilot study is presented in the next paper of this series. 
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